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1 Introduction 



The quantum theory of rotation angles (S. M. Barnett and D. T. Pegg, Phys. Rev. A, 41, 3427-3425 (1990)) is generalised to non- 

■ integer values of the orbital angular momentum. This requires the introduction of an additional parameter, the orientation of a phase 
^ , discontinuity associated with fractional values of the orbital angular momentum. We apply our formalism to the propagation of light 

^— i modes with fractional orbital angular momentum in the paraxial and non-paraxial regime. 

<n : 

> 
o 

. The orbital angular momentum of light beams is a consequence of their azimuthal phase structure [1,2]. 
Light beams with a phase factor exp(im.(/)), where m is an integer and eft is the azimuthal angle, carry 
— i ■ orbital angular momentum (OAM) of mh per photon along the beam axis [1]. These light beams can be 
^ , generated in the laboratory by optical devices which manipulate the phase of the beam, such as spiral 
phase plates or holograms. In cases where such a device generates a light beam with an integer value of 
m, the resulting phase structure has the form of |m| intertwined helices of equal phase. For integer values 
of m the chosen height of the phase step generated by the optical device is equal to the mean value of 
the OAM in the resulting beam. Recently, spiral phase steps with fractional step height [3] as well as 
£^ . special holograms [4] have been used to generate light beams with fractional OAM. Here, the generating 
| optical device imposes a phase change of exp(iM^) and M is not restricted to integer values. The phase 
structure of such beams shows a far more complex pattern; a series of optical vortices with alternating 
) i f \ charge is created in a dark line along the direction of the phase discontinuity imprinted by the optical 
£h ' device [4,5]. In order to obtain the mean value of the orbital angular momentum of these beams one has 
^ . to average over the vortex pattern. Interestingly, this mean value coincides with the phase-step only for 
q-( integer and half-integer values. Our aim is to exploit the well-known formal connections between optics 

■ and quantum theory to represent beams with fractional OAM as quantum states. We anticipate that this 
• <-h , will be particularly useful, when applied to single photons and related quantum phenomena. 

Until now, the theoretical description of light modes with fractional OAM has been based on the gener- 

■ ating optical device. For integer OAM values, however, a fundamental theoretical description exists which 
provides the only known way to treat the angle itself as quantum mechanical Hermitian operator [6]. 
This description provides the underlying theory for a secure quantum communication system [7] and it 
gives a form of the uncertainty relation for angle and angular momentum which has been confirmed in an 
optical experiment [8] . In this paper we generalise this theory to fractional values of M thereby creating a 
quantum mechanical description of fractional OAM. Such a rigorous formulation is of particular interest as 
the use of half-integer spiral phase plates has been reported recently in actual and proposed experimental 
schemes to study high-dimensional entanglement [9-11]. These studies have shown that fractional OAM 
states are characterised not only by the height of the phase step but also by the orientation of the phase 
dislocation a. For half odd-integer values of M, that is for Mmod 1 = 1/2, states with the same M but 
with a 7r difference in a are orthogonal. In light of recent applications of integer OAM in quantum key 
distribution [12] and the conversion of spin to orbital angular momentum in an optical medium [13] a 
rigorous formulation is important for possible applications of fractional OAM to quantum communication. 
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2 Fractional orbital angular momentum 



The component of the OAM in the propagation direction L z and the azimuthal rotation angle form a 
pair of conjugate variables [14]. Unlike linear position and momentum, which are both defined on an 
unbound and continuous state space, the state spaces for OAM and the rotation angle are different in 
nature; whereas the OAM eigenstates form a discrete set of states {|m)} mg z with m taking on all integer 
values, eigenstates of the angle operator are restricted to a 2ir radian interval, as it is physically impossible 
to distinguish between rotation angles differing by 2ir radians. The properties of the angle operator are 
rigorously derived in an arbitrarily large, yet finite state space of 2L + 1 dimensions [6]. This space is 
spanned by the angular momentum states \m) with m ranging from —L, —L + 1, . . . , L. Accordingly, the 
2tt radian interval [6q, 9q+2it) is spanned by 2L+1 orthogonal angle states \9 n ) with n = 9o+2Trn/(2L+l). 
Here, 9q determines the starting point of the interval and with it a particular angle operator 4>q . Only after 
physical results have been calculated within this state space is L allowed to tend to infinity, which recovers 
the result of an infinite but countable number of basis states for the OAM and a dense set of angle states 
within a 2ir radian interval. We conduct the analysis of fractional OAM in the finite dimensional space. 
After we have calculated expectation values or probabilities we can give final results in the limit of L —* oo. 

A quantum state with fractional OAM is denoted by |M), where M = m + fj, and m is the integer 
part and fx £ [0, 1) is the fractional part. We choose 'fractional' rather than 'non-integer' to follow the 
terminology in earlier work, but fi can take on every real number G [0, 1) not just rational numbers. This 
distinction, however, is not relevant for practical purposes. The state |M) is decomposed in angle states 
according to 

2L 2L 

\M) = (2L + l)- 1/2 ^exp(iM0 n )|0 n > = (2L + iy 1/2 exp(im0 n ) exp(i/i0 n )|0 n ). (1) 

ra=0 n=0 

In general exp(i/i0 ra ) will be a multivalued function but a common way to render a multi- valued function 
single- valued is the introduction of a branch cut [15], which restricts the range of the function. Usually, the 
branch cut is chosen to be along the negative real axis in the complex plane, such that the function has no 
discontinuities. Here, we choose the position of the branch cut a, that is the position of the discontinuity 
in the function exp(i/j,6 n ), independently of 9q, the starting angle of the 2ir radian interval. To achieve this 
we introduce an integer- valued function f a (9 n ) which takes on the values 1 and 0, depending on whether 
9 n is smaller or greater than a: 



fa{9 n ) '■ {9 n }n=0,...,2L ► {0, 1}, 9 r , 



1 O < On <9 + a, 
O + a < 9 n < O + 2vr. 



It is important to note that a is bounded by < a < 2tt, so that the orientation of the discontinuity is 
always understood as measured from 0o. We should stress that according to the mathematical definition 
we do not create a branch of the multivalued function, as we explicitly include the discontinuity. With this 
construction the fractional state \M) can be written as 

|M(a)> = eX ^Z^ ) ^exp(iM0 n )exp[i2vr^/ a (0 n )]|gn). (3) 
V2L+1 n=Q 

The fractional OAM state \M(a)) clearly depends on the orientation of the discontinuity a (see figure[T]), 
but the angle probability distribution is flat, P(9 n ) = (M(a)\M (a)) = (2L + l) -1 , as it is for the integer 
OAM states. For the appropriate choice of input state, phase shift and orientation this is equivalent to the 
action of a spiral phase plate operator acting on an integer OAM state [9] . It is formally possible to introduce 
a unitary operator U m ((3) which rotates the position of the discontinuity by writing U m ({3)\M(a)) = 
\M(a © (3)), where M = m = fi such that the index of the operator is matched by the integer part of 
M. We have introduced a © [3 := (a + ( 9)mod 2tt as short hand notation for the modulo 2tt addition. The 



2jtu> 



(2n-a)/d -■ 



a 



In 



Figure 1. Plot of the phase discontinuity in the fractional OAM states \M(a)). The relevant phase property is given by the factor 
cxp[ifi(8 n + 2nf a (8 n ) — a)]. For the plot we have set 0q = 0. If 9q 7^ there is an additonal global phase factor exp(i/x#o). 

action of the operator U m (P) on a general state \M'{ot)) results, not only in a rotation of the orientation, 
but also in a phase shift: 



U m {P)\M'(a)) = exp[i(m - m')(3]\M'(a © /?)). 
Further properties of this operator are described in Appendix lAl 



(4) 



2.1 Overlap of fractional OAM states 

The study of high-dimensional two-photon entanglement reported in [9-11] makes use of the fact that 
states with the same M but with a 7r radian difference in the orientation angle are orthogonal. With the 
expression for fractional states in Eq. Q we are able to calculate the overlap for two general states \M{a)) 
and \M'(a')). With help of the unitary operator U m (—a') we can set one orientation to zero and the other 
to the difference in the orientation (3 = a ( 



-a 



(M'(a')\M(a)) = {M\a')\Ul tl {-a')U ml {-a')\M(a)) = exp[i(m - m')a'](M'(Q)\M(/3)), 



2L 



( . xp [ i( „ ( m')a'] eX ^ T M ^ eX p[i(M - M')6 n ] expfi/^Tr^)], 



(5) 



2L + 



n=0 



where we have used the fact that {\6 n )}n=o,...,2L is an orthonormal set. The moduli of the rotated and 
unrotated overlaps are the same so we proceed in calculating (M'(0)|M(/3)) only. The sum over n in Eq. 
([5]) can be split in two parts by introducing an index N with On < f3 < 9n+i corresponding to the different 
cases in the definition of fp in Eq. ([2]): 



(M'(0)\M(P)) 



exp(— 
2L + 1 



N 2L 

exp(i/i27r)^exp[i(M-M')0n]+ ^ exp[i(M — M')6 n 

n=0 n =N+l 



(6) 



We find that, on substituting 6 n = 9q + 2-irn/ (L + 1) into the expression for (M'(0)\M((3)), it is possible to 
evaluate the sums using geometric progression. The overlap is a physical result and we can take the limit 



of L — > oo by expanding the exponentials in the denominators. In this limit the set of angles |#n}n=o,...,2L 
becomes dense and we can write ]im.L-+oa zZ'rp' = nm L^oo &N+1 = 

(M'(0)|M(/3)) = exp(-i^) ieX 2 P ^~^^ o] (exp^vr) (l - exp[i(M - M')(3]) 
+ (exp[i(M - M')/3] - exp[i(M - M')2tt]) } . 

From this expression it is possible to calculate the modulus square of the overlap {M'(0)\M(fl)) with help 
of several trigonometric identities: 

\(M'(0)\M(P))\ 2 = {M _ l MI) 2^ [^(Mvr) +sin 2 (M'7r) 

- 2cos[(M - M')(n - (3)} sin(Mvr) sin(M'7r)] . 

This overlap probability is plotted for different values of (3 in figure [2j For (3 = the plot shows that 
there are parallel lines of equal overlap probability for the diagonals with M — M' = K £ M. The overlap 
probability along these lines is given by sin 2 (tt K) / (it 2 K 2 ) . For j3 ^ the overlap probability along the 
diagonal lines M — M 1 = K becomes dependent on M and M' . In particular for (3 = n and M = M' the 
overlap probability reaches zero for fj, = jj! = 1/2. This was first reported in [9] and, if we specialise the 
expression for the overlap probability in Eq. j7| to M = M' , we obtain an identical expression for the 
overlap probability {M (0)\M' '{(3)) which only depends on the fractional part \i and the orientation (3: 

\{M(Q)\M{(3))\ 2 = (\-^sin 2 (^7r) + cos 2 (/i7r). (9) 

This property has been used in an experimental detection of entanglement for fractional orbital angular 
momentum [11] and in a proposed experiment to test for non-locality of orbital angular momentum states 
[10]. 



2.2 Angular momentum distribution for fractional states 

We can decompose a fractional OAM state into the basis of integer OAM states. The amplitudes of the 
decomposition are given by the expression for the general overlap ([7|), where we choose M' = m' G Z, and 
we obtain: 

c m ,[M(f3)\ = <m'|M(/?)> = expC-i^) 16 ^ 1 ^"^^ [exp[i(m - m')(3] (1 - exp(i/,2vr))] . (10) 

Only the complex argument of the probability amplitudes c m >[M((3)] depends on the relative orientation; 
the probabilities P m >(M), given by the modulus square of the probability amplitudes, are independent of 
0: 

Integer OAM states form an orthonormal basis, that is they obey (m\ 77i/) — $mm' fo^" Tfl^ Tflf G Z. This 
particular result is recovered for the OAM probabilities P m >(M) from Eq. (I10p . For M = m £ Z and 
m m' the fractional part fi is zero resulting in a vanishing OAM probability. For M = m', we can 






Figure 2. Plot of the overlap probability for two general OAM states \M) and \M') for three different values of the difference in the 
orientation f3. a) /3 = 0: Along parallel diagonal lines, for M — M' = K G Z the overlap probability is independent of the particular 
values M and M' . b) fi = 7r/2: The overlap probability generally depends on the fractional values, but does not reach zero for M = M' . 
c) /3 = 7r: For M = M' and fj, = /j,' = 1/2, the overlap probability reaches zero. 



determine the value of P m >(m ) by the limiting procedure: 



P , (m) = lim P m , (m + fi) = — ^ lim — COb ( 2 / i7 ^) = $ _ ( 12 ) 
m^o 27T 2 /i->o (to + // — m/) 2 

For integer values of M the OAM distribution is thus singular, consisting of a single non-vanishing prob- 
ability at M = m' . For fractional values of M, however, the probabilities are peaked around the nearest 
integer to M, as can be seen in figure O 

Owing to the completeness of the OAM basis states \m'), the probabilities P m >(M) sum to unity. The 
relevant summation can be executed using the contour integration method [16]. In a similar way the mean 
value of the OAM can be calculated: 

v 2 ^ /„ ,wn l-cos(2Mvr) ^ m' _ sin(2Mvr) . . 

m'=— co m'=— oo v ' 

where we have used cos(2^7r) = cos(2M7r). The result is an OAM mean that is equal to M only at integer 
and half- integer values of M (see figure 13} • This expression was first reported by Berry in [4] where it was 
also confirmed by experiments. The OAM variance, however, tends to infinity. From a more general point 
of view it is not surprising that the OAM variance is divergent. The discontinuity in the phase leads also to 
a discontinuity in the angle wavefunction. For such a discontinuous wavefunction the mean square OAM 
will be divergent [17]. 



3 Propagation of a fractional beam 

The decomposition of a fractional state in OAM eigenstates can be applied directly to the decomposition 
of a light mode emerging from a fractional phase step into OAM eigenmodes. By using a decomposition 
into Bessel beams rather than Laguerre-Gaussian beams we can switch conveniently from a non-paraxial 
solution to a paraxial solution by changing the propagation factor in the Bessel beam. Berry has reported 
the propagation of a plane wave with a fractional OAM imprinted on it [5]. Here, we find many of the 
same features for a beam of finite transverse extent. 

For the decomposition we consider a Gaussian input beam with a finite beam width wq. The beam 
emerging from spiral phase plate can thus be written as 

*(p, <p,z = 0) = exp (-^i ) exp(iM<^). (14) 

A further advantage of Bessel functions as solutions to the exact and paraxial wave equation is that the 
dependence on the radial and longitudinal coordinates factorises. We can therefore express the Gaussian 
beam directly behind the phase step as a superposition of Bessel beams with different transverse wavenum- 
bers k. The propagated wave is obtained by multiplying each term in the superposition by the appropriate 
propagation factor for the exact or paraxial theory. The superposition of Bessel beams is specific for every 
integer value of the orbital angular momentum and the decomposition of the wave behind a fractional 
phase step requires a summation over integer values m weighted by the coefficients c m '[M(/3)]: 

oo „ 
*(p, ip, z = 0) = ^2 c m ,[M(P)]exp(im'tp) dKd m/ (K)J m ,(Kp). (15) 

m'=— oo 

The coefficients c m > [M(/3)] are given by the overlap between the integer and fractional OAM states 
{m'\M{(3)) (see Eq. fjlOf) ) . This follows on multiplying both expressions for ty(p,(p,z = 0) in Eqs. (I14p 
and (|15p by exp(— im"4>) and integration over ip. The radial part contains the function d m i{K) which is 
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Figure 3. Plot of the OAM distribution of fractional states for different fractional values fi. The distribution has a peak at the nearest 

integer to M. Choosing a different integer value m! shifts the distribution uniformly by m' — m, such that the distribution is peaked 
around m' instead of m. In b), for fi = 0.5 this results in two peaks of equal height at the two neighbouring integers. The spread in the 

distribution is determined by the fractional value fj,. 
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Figure 4. Plot of the OAM mean value for fractional OAM states. Only at integer and half integer values of M is the OAM mean M 

equal to M. In general the relation between M and M is given in Eq. (1 1 3 D - 



determined by the integral 



p 2 



6XP V 2u^J = J dKd m'(^)Jm'(Kp). (16) 
On using the Fourier-Bessel theorem [18], the function d m t(K) can be found to be: 

P 2 
'o 



d m >(K) = n I dp pJ m >(Kp)exp ( -7^2 ) • ( 17 ) 



The solution of this integral can be written in terms of modified Bessel functions I u with a half-integer 
order v [19]. For all positive half- integer orders the modified Bessel functions increase exponentially. But 
this behaviour changes for v = — ^ and v = | which renders the integral expression for the coefficients in 
terms of modified Bessel functions invalid for m! < —2. We therefore have to evaluate the integral in Eq. 
(fT7|l for negative values of w! separately by using J_ m '(Kp) = (— l) m J m >(Kp) [16]. The coefficients in Eq. 
(I17p may thus be written as 

(_l)Kl m '<0, 

m'>0. { ' 



- (***«p (-T 4 ) (-T 4 ) " ^ (-T 4 )] * { 1 



On substituting the coefficients in Eq. (|15j) the wave emerging from a fractional phase step is determined 
at z = 0. To calculate the propagated wave we have to add the appropriate propagation factor in the integral 
in Eq. (|15|) . For the exact form the propagation factor is given by exp(i\/& 2 — k 2 z) and the propagated 
wave reads as 

^(p,ip,z)= ^2 c m <[M(/3)] exp(imV) / d«; d\ m ,\ («) J\ m ,\ exp(i\/ A; 2 - k 2 2;). (19) 

m'=-oo ^° 

A distinction between negative and positive values of m' is not necessary in this formulation as for negative 
m! the alternating factor (— l) m is compensated by the alternating sign of the Bessel functions J- m ' ( k p) = 



(— l) m J m >{Kp). The integration boundaries are set from zero to infinity and this includes evanescent waves 
for k > k. As pointed out in [20], the evanescent components are necessary to describe the singularity 
correctly. It is worth noting, however, that the character of the integral changes at k = k, where the square 
root turns imaginary and the exponential changes from an oscillating to a decreasing behaviour. In the 
paraxial solution, where the propagation factor is exp[ifc(l — l/(2n 2 )z], this distinction does not exist. The 
z dependent term exp(ifoz) in the propagation factor can be written in front of the integral, which gives 
for the paraxial solution 

ip(p,<p,z)= 22 Cm'[M(P)]exp(im'(p)exp(ikz) d« d| m /|(«) J\ m ,\{Kp) exp ( -i-^z J . (20) 

m'=-oo ^° \ K / 

The integrals for the exact and paraxial solution are solved numerically to give the intensity and the phase 
profile of the propagated wave. 



4 Numerical results 

For the numerical calculation we only take eigenmodes into account for which the corresponding OAM 
probability |c m '[M(/3)]| 2 is larger than 10 -4 . The widest spread in the OAM spectrum occurs for half odd- 
integer values of M and for these the restriction to |c m <[-^(/3)]| 2 > 10 -4 amounts to 66 integer eigenmodes 
contributing to the decomposition. The phase and intensity profiles are calculated for the fractional phase 
step M = 3.5, because for half odd- integer phase steps we expect to find the interesting phenomenon of 
the formation of a chain of vortices in the line of low intensity created by the phase discontinuity [4,5]. We 
give the results for both the exact solution and in the paraxial approximation. We choose a beam width 
of wo = 10000 \/(2ir), which corresponds to a beamwaist of about one millimeter for visible light. For 
the paraxial solution and the chosen beam width the intensity profile remains unchanged on propagation 
and the phase profile simply rotates around the propagation axis. This is why we present four profiles 
for the exact solution at the propagation distances kz = 1,5,50,200 and two for the paraxial solution at 
kz = 1, 200. The results are presented in two series of graphs in figures [S] and In the intensity profiles the 
most obvious feature is the radial line of zero intensity. For a value of M with a fractional part different 
from 1/2 the intensity will be low but not equal to zero along this line. The other important feature in 
the intensity profiles is the number and position of the spots of zero intensity in the centre of the beam. 
The graph for M = 3.5 shows 3 such spots. These spots correspond to optical vortices and the number 
of vortices is given by the modulus of the nearest integer to M. While a beam with integer OAM m 6 Z 
propagates with an optical vortex of charge m on the axis, beams with fractional OAM only show vortices 
with charge ±1. None of these vortices is on the axis, but the whole central region has low intensity. 
Another prominent feature are the diffraction fringes surrounding the axis; these arise due to diffraction 
of the initial singularity at the centre. The radial fringes in the direction of a are due to diffraction on the 
discontinuity in the spiral phase plate or hologram. The contrast in the fringes is highest for a fractional 
part p = 1/2. It is diffraction, of course, which is responsible for breaking the initial rotational symmetry. 

The phase profiles show clearly that the fractional strength vortex splits up in vortices with strength ±1. 
The phase is represented in false colour, such that the colours for the phase of and 2ir are identical. A 
vortex corresponds to a point where all colours meet and for a vortex of strength ±1 a circle surrounding 
the vortex goes through the colour circle once. The sign of =Fl determines if the colours run through 
the circle clockwise or anti-clockwise. The most prominent feature in the phase profiles is the creation of 
additional vortices in the region of zero intensity for half odd-integer phase steps. In figure ([6k), for kz = 1, 
one can see that there are no additional vortices, while in figure ([6h>) one pair of vortices has formed and 
another pair is in the process of forming. At greater propagating distances the series of vortex pairs moves 
closer together. From the different rotation sense for the colours it can be seen that the vortices in this 
series have alternating charge. The vortices annihilate each other differently depending on whether the 
phase step has a fractional part p of just under 1/2 or just over 1/2 [5]. This leads to the formation of 
a new vortex or the destruction of an existing vortex for fractional phase steps with p > 1/2. The phase 



Exact solution, M=3.5, k?=\ Exact solution, W=35, kz=5 




-111 -S -o -4 -2 2 -I f P 10 -20 -15 -10 -5 5 10 15 20 



Exact solution. M=3S, k?.=50 Exact solution. M=35. k/=2P0 




-to jo -20 o 2fl 40 60 -ion -?o a so im 



Paraxial solution, M=3.5, kz^l Paraxial solution, M=3.5, k>=200 




-lO-S-ci-4 -2 2 -I 6 P 10 -ID0 -50 50 100 



Figure 5. Intensity profiles for M = 3.5. The graphs a-d) show the exact solution for the propagation distance kz = 1, 5, 50, 200 
respectively. The graphs e) and f) show the paraxial solution for the distances kz = 1 and kz = 200. For fractional phase steps with 

/i = 1/2 the intensity along the radial line drops to zero. 
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Figure 6. Phase profiles for M = 3.5. The graphs a-d) show the exact solution for the propagation distance kz = 1, 5, 50, 200 
respectively. The graphs e) and f) show the paraxial solution for the distances kz = 1 and kz = 200. The phase profiles show the 

creation of additional vortices on propagation. 



pattern rotates on propagation. For a light beam with integer OAM this rotation would be uniform, but 
due to the presence of the imprinted phase discontinuity, which does not change under propagation, the 
rotation of the phase pattern changes. If one follows for example the red line of zero phase roughly along 
the negative y-axis in figure [6^, this line will rotate counter-clockwise on propagation and at some distance 
it will be oriented along the negative x axis. Within the same propagation distance the line of zero phase, 
that is roughly oriented along the positive y axis, will rotate through the phase discontinuity until it is 
oriented along the negative y axis. Over the same propagation distance these two lines of zero phase cover 
different rotation angles. This shows that the phase discontinuity at the angle a causes the non-uniform 
rotation of the phase pattern. 

On propagation fractional modes develop a number of features which can only be seen in the exact 
solution. Although the paraxial solution approaches the exact solution for z — > oo the paraxial solution 
does not show the changing position of the integer vortices or the formation and movement in the chain 
of vortices. 



5 Conclusion 

In this article we have derived a rigorous quantum formulation of fractional orbital angular momentum 
(OAM). This formalism is a generalisation of the quantum theory of rotation angles by Barnett and Pegg [6] 
to fractional values of the OAM. Fractional values of OAM require the introduction of a branch cut in 
the angle representation of the OAM eigenfunctions. The orientation of the branch cut is an additional 
parameter in the description of fractional OAM. We have calculated the overlap of two general fractional 
OAM states and we have found that for half odd-integer values of the OAM and a relative orientation of 
7r radian the overlap is zero and the corresponding states are orthogonal. This confirms earlier work on 
fractional OAM with half- integer spiral phase plates [9]. 

We have applied our theory to the propagation of light beams with fractional OAM. For light carrying 
integer OAM a stable optical vortex of corresponding integer strength forms on the propagation axes. For 
fractional OAM no fractional strength vortices are stable on propagation, and instead a number of strength 
±1 vortices are formed in the central region of the beam. Additionally, a line of low intensity is visible 
corresponding to the orientation of the branch cut or the phase step discontinuity in the generating optical 
device. For half-integer phase steps a chain of additional vortices with alternating charge is formed in this 
region of low intensity. We have confirmed theoretical and experimental results regarding the evolution of 
optical vortices in light beams with fractional OAM. Further to existing theories we have calculated the 
propagated fields in the paraxial and non-paraxial regime for an incident beam of finite width. The gradual 
formation of additional vortices in the region of low intensity is a phenomenon which is only visible in the 
non-paraxial solution. 
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Appendix A: Rotation of the orientation of the discontinuity 

Quantum states with fractional OAM depend on a, the orientation of the phase discontinuity. To rotate 
the orientation we introduce a class of operators U m (f3),m 6 K,/j 6 [0, 2tt) defined by their action on a 
state with fractional OAM: 



U m {p)\M'{a)) = exp[i(m - m')0\\M(a /?)), 



(Al) 



where (3 is the angle through which the discontinuity is rotated. The addition a © p = (a + pjmod 2ir 
yields a result in the range [0, 2tt). For m = m', that is if m is equal to the integer part of M, this operator 
is a pure rotation of the orientation a, but for m ^ m' the rotated state acquires also a phase shift unless 
/3mod2-7r = 0. Consequently, the operator acts as identity operator if (3 is an integer multiple of 2tt. The 
eigenstates of this operator are the integer OAM states. Acting with U m {(3) on an integer state \m),m G Z, 
results only in a phase shift: 

U m (/3)\m') = exp[i(m - m')(3]\m) . (A2) 
All the eigenvalues have unit modulus and the operator U m {(3) is therefore unitary with U m ((3) = 

u m {-(3) = ul{p). 

From the definition of the operator in Eq. (lAlj) it is obvious that the set of operators cannot form 
a group under multiplication for arbitrary combinations of parameters (3 and m. But for a fixed m the 
set of operators {Um(p')} 0e[o,2w) does form a one-parameter group under multiplication with the group 
parameter (3 [21]. As proof we show that U m ((3')U m {(3) can always be written as Umil) with 7 6 [0, 2tt): 

U m (f3')U m {(3)\M'(a)) = U m {(3>) exp[i(m - m')(3]\M'(a © (3)), 

(A3) 

= exp[i(m - m')(J3 + fi)]\M '{a ®(3® /?')) = U m {(3 © /?')|M'(a)>- 

Therefore U m {(3')U m {(3) = U m {(3 © /?') as the modulo 2ir addition (3 © is always in the interval [0, 27r). 
There is also a neutral element in form of the identity operator for (3 = 0. Moreover, to every operator 
U m ((3) exists an inverse element in the form U m {(3) = U m (—(3). 

The product of two operators with the same orientation but different values for the fractional orbital 
angular momentum can be combined to give 

U m ((3)Um'((3) = exp[-i(m + m')(3}U m+ml ((3 © (3). (A4) 

The set of operators {C/ m (/3)} m ez does not form a group under multiplication for (3 ^ 0. For (3 = these 
operators are all identity operators. 

We have introduced the operator U m (/3) to facilitate the calculation of the general overlap (M' (a')\M(a)} 
in section I2TT1 The mathematical properties suggest, that the action of this operator on a state \M'(a)} 
is independent of the fractional part //. A physical implementation of this operator would be an optical 
device with two edge dislocations; one to compensate the original dislocation at an angle a, the other to 
imprint the new phase discontinuity at an angle a (3. This implementation, however, would be specific 
for given fractional part // . As the effect of the edge dislocation is completely specified by the fractional 
part fi' and the orientation a it would be possible to change the orientation of the state \M'(a)) with an 
implementation of the operator U m ((3), where the integer part of M' is different from m, which causes the 
overall phase factor in Eq. (|Aip . 
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